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Conscan is commonly used by DSN antennas to allow adaptive fracking of a 
target whose position is not precisely known. This article describes an algorithm 
that is based on a Kalman filter and is proposed to replace the existing fast Fourier 
transform --based (F FT- based) algorithm for conscan. Advantages of this algorithm 
include better pointing accuracy , continuous update information, and accommoda- 
tion of missing data. Additionally, a strategy for adaptive selection of the conscan 
radius is proposed. The performance of the algorithm is illustrated through com- 
puter simulations and compared to the FFT algorithm. The results show that the 
Kalman filter algorithm is consistent]}' superior. 


I. Introduction 

The objective of conscan is to improve the estimate of 
the target location; the target location is defined as the 
direction of maximum carrier power as received on the 
ground. During conscan, the antenna is rotated in a circle 
about a point we call the conscan center at a constant rate 
w. Figure 1 shows the variables used in the analysis. 

Without loss of generality the origin is placed at the 
conscan center. We define x = x(t) as the 2- vector which 
specifies what we seek to estimate: the unknown actual 
target at time t y and we define x = x(/) as the estimate 
of the target’s location. We assume that discrete mea- 
surements are taken at times t — Let the subscript “z” 
denote a quantity taken at time ti \ for example, x, = x(fj). 
Lack of the “i” subscript implies that the affected quantity 


is constant over a conscan period. We refer to x; as the 
state of the system, and thus x* is the state estimate. The 
known instantaneous location of the antenna boresight is 
defined by x al} which during conscan is given by 


x 


ai — 


R cos (a?/,;) 
R sin (udj) 


(i) 


Here, R is referred to as the conscan radius, and w is 
the conscan frequency. 


The output of the antenna receiver is the carrier power 
Pqi . Assuming a circularly symmetric antenna gain func- 
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tion , the carrier power can be approximated as a quadratic 
function of the offset angle: 1 

P ci = P oi (l - + T}i (2) 

where P 0 i is the maximum carrier power which occurs 
when the antenna is pointed directly at the target, h is 
the antenna half-power beamwidth, /i = 4 In (2), /?,* is the 
offset angle between the antenna boresight and the target, 
and rji is the signal noise, or error. We assume that the 
signal noise is random with zero mean and is Gaussian dis- 
tributed with a standard deviation of a v . From Fig. 1, we 
have 


Pi = - x ai ) T (xi - x a! ) (3) 

and thus 


and thus 

2P T 

Zi = -^-x at Xi + rti ( 8 ) 

In the subsequent analysis, we refer to z t as our mea- 
surement. 

II. FFT Estimator 

A brief discussion of the existing fast Fourier transform 
(FFT) estimator is given below so that it can be compared 
to the Kalman filter estimator. We assume that Po* = Po 
and x, — x are constant over a conscan period. Multiply- 
ing Eq. (8) by x a * and averaging both sides over a conscan 
period yields, after some algebra, 

( ZiXai ) ~ X + (rjiXai) (9) 


P«(x,) - Po.' (l - - x al ) T (x, - x a! )) + i n (4) 


Utilizing the fact that rji is assumed to be zero mean, 
we obtain an estimate for the state: 


If we assume that P 0t - = Po and x, = x are constant 
throughout a conscan period, then from Eq. (1) it is easily 
shown that the ensemble average of the carrier power over 
the conscan period is approximately 


h 2 (Z|X QI ) 

PofiR 2 


( 10 ) 


We define M to be the covariance of the state estimate: 


(P c> ) « P 0 ( 1 - ^(x T x + R 2 )) ( 5 ) 


M sP((x i -x)(x,-x) T ) (11) 


Subtracting Eq. (5) from Eq. (4) yields 


where E denotes expectation, or ensemble average. After 
further manipulation it can be shown that 


Pci(Xi) - (Pci) » ^p^xJ;Xi + T)i (6) 

Equation (6) constitutes the measurement equation for 
the purposes of the estimators. It is important to note 
that the measurement equation is linear with respect to 
the unknown state x* } as this makes it possible to apply 
linear estimation theory to estimate x t . 

We define z * as 

Zi = P ei (xi) - (Pci) (7) 


3 L. Alvarez, “Open Loop Conscan Pointing Error Estimation Accu- 
racy at Ka-band ” JPL Interoffice Memorandum 3328-93-044 (in- 
ternal document), Jet Propulsion Laboratory, Pasadena, Califor- 
nia, July 8, 1993. 


M = 


2Pq fi 2 R 2 n 


I 


( 12 ) 


where I is the 2x2 identity matrix and n is the number 
of samples taken over a single conscan period. This result 
is consistent with results described elsewhere. 2 

In practice, an FFT is used to estimate x [1]; this yields 
the same result as Eq. (10) but also yields harmonics above 
the fundamental frequency. Hence, we call this the FFT 
algorithm. Use of an FFT imposes several constraints on 
the algorithm: n must be an integer power of 2; .there 
must be no missing measurements; and the measurements 
must be evenly spaced in time. Theoretically, failure in 


2 Ibid. 
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the estimator can be detected by the presence of signif- 
icant amplitudes in the harmonics above the fundamen- 
tal frequency, but in practice, detection of the harmonics 
above the fundamental frequency has been only partially 
implemented, partly because it is readily apparent when 
conscan has failed. 


III. Kalman Filter Estimator 

A detailed general description of a Kalman filter is 
given in [2]. However, square root algorithms in which 
the Cholesky factors of the covariance matrix are prop- 
agated are widely recognized to have superior numerical 
properties, and hence we implement the square root co- 
variance filter described in [3]. We assume that between 
time intervals ti and x; drifts according to 


Xi+i=Xi+ni (13) 


where n* is a random 2-vector with covariance Q = 
E(uinJ). We assume that at time t{ the state estimate 
is known with covariance M,- ~ E ((x, — x*)(xi — x;) T ). 
Our measurements take the form 


z i = H iXi + v* 


(14) 


<*>,= ; t * 


fc=ri — n + 1 


Poi = 




1 - -f^(R 2 + x, T x.) 


Z, = Pci-(Pc). 


Xai — 


R cos (u )t{) 
[Psin (ujti) 


Hi = 


T = 


2PoiA* T 


T x U i = 


R 1 / 2 HiS, 0 
0 Si Q 1 ' 2 


v t 1/2 0 0 

L Gi S; + 1 0 


K, = G,v: 
M , = S,S, T 
Zi = H ,x, 


1/2 


(16) 


where v* is a random vector with zero mean and whose 
covariance is given by R = E ( v < v 7). From Eq. (8), we 
can make the identifications 


H, = 


2Pp/i T 

h 2 ai 


Vi = T]i 


R= a 2 ,, 


(15) 


With these preliminaries completed, we are ready to 
construct the Kalman filter estimator, which requires ex- 
ecution of the following updates at each time interval 


x,*+i - x t - + K,(zi - Zi) J 

Here, n is the number of samples in a complete con- 
scan period. Uj is an orthogonal transformation that tri- 
angularizes T and may be obtained, for example, using 
Householder transformations. Quantities such as R 1 / 2 are 
“matrix square roots” or Cholesky factors. The Kalman 
filter algorithm is straightforward except for the estimation 
of the average carrier power, which is computed as 

a running average of the n most recent outputs from the 
antenna receiver, and estimation of the maximum carrier 
power, Poi> which is computed using Eq. (5). Thus, the 
Kalman filter requires that at least n measurements have 
been accumulated before the filter can be started. This 
results in a delay of a single conscan period before the 
Kalman filter begins producing estimates. A similar delay 
also occurs in the FFT algorithm. 
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IV. Selection of the Conscan Radius 

An algorithm for selecting the conscan radius was de- 
veloped as part of this study. We again assume a constant 
target offset x, = x. Recall from Eq. (5) that the average 
carrier power over a conscan period is given by 

<P c ) = Po(l-£(x T x + / 1 ! 2 )) (17) 

Obviously, the average carrier power is maximized for 
x = 0, which can be attained by repointing the antenna 
until the conscan center coincides with the target loca- 
tion. However, the best that can be done in practice is 
to repoint the antenna to place the conscan center at the 
estimated target location x. In the newly pointed antenna, 
the average carrier power becomes 

(Pc) = Pa (l - ~ ((x - x) T (x - x) + fl 2 )) (18) 

and the average-carrier-power expectation, using the defi- 
nition of M given in Eq. (11), becomes 

E({P C )) = P 0 (l - ~ (trace(M) + 7? 2 )) (19) 

A necessary condition for the average power to be max- 
imized is for the expectation of the carrier power given by 
Eq. (19) to be stationary with respect to R ; taking a dif- 
ferential and setting it to zero requires 

— trace(M) = -2 R (20) 

An analytic expression for M (R) for the Kalman filter 
would be difficult if not impossible to obtain. In princi- 
ple it would be possible to propagate d^A/dR along with 
the Kalman filter equations, but this more than doubles 
the number of computations performed at each iteration 
and excessively complicates the algorithm. Instead we use 
Eq. (12), developed from FFT principles, as an approxima- 
tion for M, with the reasoning that this result is an upper 
bound on M because it is based on data from only a single 
conscan period. In this case, solving Eq. (20) yields the 
optimal conscan radius: 



Typically the conscan radius is set to result in 0.1 dB 
gain loss in the carrier power. The set of parameters sum- 
marized in Table 1 yields an optimal conscan radius of 
R op t — 1 -86 mdeg for a 34-m antenna. The conscan radius 
for 0.1-dB gain loss is 5.9 mdeg. This suggests that the 
conscan radius is currently larger than it needs to be, at 
least for this set of parameters. 


V. Implementation Details 

The behavior of the Kalman filter depends heavily on 
the selection of the process noise and measurement noise 
covariances and also on the covariance of the state estimate 
used to initiate the algorithm. The following criteria are 
used to select these values. 

The measurement noise variance R is assumed to 
be provided by the RF receiver. (See [4,5] and other 
material 3,4 for information.) For the simulations, a model 
developed elsewhere 5 is used. Though R is not indepen- 
dent of offset angle, it can be considered to be constant for 
the small offset angles used during conscan. 

The process noise covariance Q can be established by 
considering the maximum anticipated drift in the offsets 
within one sample period. We select 

Q = (drift rate x sample time) 2 I (22) 

where I is the 2x2 identity matrix. The drift rate used 
in the simulations is 

, . r . x i r 

drift rate = - x : — (23) 

5 conscan period v 

An alternative, which we suggest for implementation, 
would be to use the reciprocal of the antenna gimbal servo 
bandwidth. 

The initial covariance M 0 of the state estimate is given 
hy the expected uncertainty of the initial target location. 
This should be based on the known open loop (blind) 

3 M. Aung and S. Stephens, “Statistics of the P c /N 0 Estimator in 
the Block V Receiver,” JPL Interoffice Memorandum 3338-92-089 
(internal document), Jet Propulsion Laboratory, Pasadena, Cali- 
fornia, April 29, 1992. 

R. E. Scheid, “Statistical Analysis of the Antenna Carrier Power,” 
JPL Interoffice Memorandum 343-92-1291 (internal document), Jet 
Propulsion Laboratory, Pasadena, California, October 9, 1992. 

5 Ibid. 
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pointing performance of the antenna. For the simulations, 
we assume that the target can be located anywhere within 
the conscan circle; thus we choose 

M 0 = R 2 I (24) 


VI. Computer Simulations 

The parameters used in the computer simulations are 
typical for a 34-m DSN antenna operating at Ka-band 
(32 GHz) and are summarized in Table 1. A random- 
number generator was used to simulate noise in the carrier 
power. Each simulation used a different seed to start the 
random-number generator, which ensures that a reason- 
able cross section of cases is represented. The Kalman 
filter estimator was then propagated through the data, 
starting with a zero initial estimate. For comparison to 
the existing algorithm, the FFT estimator was also prop- 
agated. The FFT estimator is easily distinguished from 
the Kalman filter because updates are produced only once 
per conscan period. The simulations were written in Mat- 
lab, which allows symbolic manipulation of matrices and 
vectors within an interpretive programming environment. 

The first simulation is shown in Figs. 2- 4 and illustrates 
the convergence of both Kalman filter and FFT estimators 
to a constant target offset. The simulated carrier power is 
shown in Fig. 2. The sinusoidal modulation of the carrier 
power resulting from the nonzero offset is readily appar- 
ent. Figure 3 shows the output of the Kalman filter and 
the FFT estimator. The dashed line on the plot indicates 
the true target location. The Kalman filter estimator be- 
gins by collecting carrier pow r er over a single conscan pe- 
riod so that an estimate of the average carrier power (P C i) 
can be obtained. Propagation of the Kalman filter begins 
immediately thereafter. Convergence of the estimate is 
virtually complete within a fraction of a conscan period. 
The Kalman filter estimates are consistently better than 
the FFT estimates. In fact, in steady state the standard 
deviation of the Kalman filter estimates is about one-third 
that of the FFT estimates. 

A measure of the convergence of the Kalman filter is 
given by the norm of the estimate covariance M*. Fig- 
ure 4 shows the time history of |M,| and demonstrates 
that the Kalman filter converges very rapidly at first and 
then reaches steady state after about 100 sec, or about 
3 conscan periods. This is an important observation be- 
cause it assures us that running the simulations longer will 
not reveal any different behavior of the algorithm. 


In the simulation shown in Fig. 5, the measurement 
noise has been doubled to investigate the relative sensitiv- 
ities of the two algorithms to measurement noise. Con- 
vergence of the Kalman filter algorithm is changed very 
little, though close examination reveals a larger amplitude 
in the high-frequency component of the Kalman filter esti- 
mates. Errors in the FFT estimates are pronounced, being 
roughly doubled when the measurement noise is doubled. 
As before, the Kalman filter estimates are about three 
times more accurate than the FFT estimates, as measured 
by the standard deviations in steady state. 

The simulation shown in Fig. 6 shows the ability of 
both algorithms to successfully estimate a relatively large 
initial offset. Despite the fact that the initial offset was 
set at more than three times the conscan radius, both al- 
gorithms quickly converge to the correct value. In general, 
convergence is improved in both algorithms as R grows 
larger. 

In Fig. 7, a data dropout between 100 and 150 sec was 
simulated. The measurement covariance update and state 
update steps in the Kalman filter algorithm are not applied 
for each missing measurement. The most recent estimate 
is used throughout the dropout period, and then when 
data become available again, the Kalman filter gracefully 
continues. Recall that this FFT algorithm is incapable of 
accommodating missing data, so there is no corresponding 
plot showing FFT estimates. 

The ability of the algorithms to accommodate a drifting 
parameter is illustrated in Fig. 8. For this simulation, the 
offset is taken to be a constant drift at a rate of 2 mdeg 
per 100 sec. Both algorithms successfully track the offset, 
though they consistently underestimate the offset slightly. 

The performance of the algorithms in the presence of 
a gradually changing maximum carrier power P 0 is shown 
in Fig. 9. In this simulation, the nominal carrier power 
is increased by 2 percent over each conscan period. The 
Kalman filter algorithm behaves roughly as before, though 
some oscillation is apparent in the estimates. Numerical 
experiments have shown that the amplitude of these oscil- 
lations increases with the drift rate Pq, but degradation of 
the Kalman filter estimates is gradual and, in this case, is 
comparable to the FFT algorithm estimates. 

VII. Summary, Conclusions, and Future 
Work 

A sequential estimator for executing conscan on the 
DSN antennas has been described. The algorithm, based 
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on a Kalman filter, has important advantages over the 
existing FFT algorithm, including better accuracy of the 
estimates, nearly continuous updates, and the ability to ac- 
commodate missing data. It also has many of the desirable 
properties of the FFT algorithm, particularly robustness 
and the ability to track large initial offsets. In addition, 
an algorithm that is used for selecting the optimal con- 
scan radius and is based on maximizing carrier power is 
proposed. The performance and behavior of the Kalman 


filter algorithm as compared to the FFT algorithm have 
been analyzed extensively through computer simulations, 
which have shown that the Kalman filter algorithm is con- 
sistently superior. 

Future work on the Kalman filter algorithm should con- 
centrate on a field implementation to demonstrate that it 
can work as well with real data as it works in an idealized 
computer simulation. 


References 


[1] J. E. Ohlson and M. S. Reid, Conical-Scan Tracking with the 6$-m Diameter 
Antenna at Goldstone } Technical Report 32-1605, Jet Propulsion Laboratory, 
Pasadena, California, October 1, 1976. 

[2] A. E. Bryson and Y. Ho, Applied Optimal Control , New York: Hemisphere Pub- 
lishing Corporation, 1975. 

[3] M. Verhaegen and P. Van Doorcn, “Numerical Aspects of Different Kalman Filter 
Implementations,” IEEE Transactions on Automation and Control , vol. AC-31, 
no. 10, pp. 907-917, October 1986. 

[4] T. J. Brunzie, “The Parkes Front-End Controller and Noise-Adding Radiome- 
ter,” The Telecommunications and Data Acquisition Progress Report $2-102, 
vol. April- June 1990 , Jet Propulsion Laboratory, Pasadena, California, pp. 1 19— 
137, August 15, 1990. 

[5] A. M. Monk, “Carrier- to-Noise Power Estimation for the Block V Receiver,” The 
Telecommunications and Data Acquisition Progress Report $2-106, vol April- 
June 1991 , Jet Propulsion Laboratory, Pasadena, California, pp. 353-363, Au- 
gust 15, 1991. 


228 



Table 1 . Selected parameters used in the simulations. 


Quantity Value 


Description 


n 

t 

R 

R 

M 0 

Q 


h 

Po 


32 

1 sec 

5.3 X 10“ 15 W 


5.9 mdeg 

R 2 1 



Number of samples in a conscan period 

Sample period 

Standard deviation of Pq 

Measurement covariance 

Conscan radius for 0.1-dB gain loss 

Starting state covariance 

Process noise covariance 

Starting estimate 


65 mdeg 
4.14 x 10“ 13 W 


Half-power beamwidth 
Maximum carrier power 


y 


t 



Fig. 1. Definition of variables. 


ESTIMATED ELEVATION OFFSET, mdeg | CARRIER POWER, Wx 10~13 



TIME, sec 

Fig. 2. Simulated carrier power used for testing conscan 
algorithms. 



Fig. 4. Convergence of the estimate covariance for a constant 
offset. 




Fig. 3. Convergence of Kalman filter and FFT estimates to a 
constant offset. 


Fig. 5. Performance of the estimators for doubled measurement 

noise. 


230 





ESTIMATED ELEVATION OFFSET, mdeg ESTIMATED ELEVATION OFFSET, mdeg 



Fig. 6. Tracking of a large offset. 



Fig. 8. Tracking a constant rate drift in the offset. 



Fig. 7. Demonstration of the Kalman filter algorithm with a data 
dropout. 



TIME, sec 

Fig. 9. Effect of maximum carrier power drift on estimator per- 
formance. 
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